Collective Modes and Emergent SO (6) Symmetry in the Iron Pnictides 
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We show the existence of an emergent S0(6) symmetry in the low energy description of the iron 
pnictides. This approximate symmetry provides a unifying framework for the occurrence of spin 
density wave (SDW) and superconductivity (SC) in these materials. We use this symmetry to make 
several predictions for future experiments, including the topology of the phase diagram and the 
presence of various resonant modes in neutron scattering experiments in both the SC and SDW 
phases. We also predict the existence of a new "Orbital Density Wave" state, which competes with 
both SDW and SC orders. 
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Introduction. — The recent discovery of superconduc- 
tivity in LaOFeP [1] and subsequently in LaOi_a;Fj.FeAs 
[2] has generated great excitement in the iron pnictides, 
a new class of materials with great potential for applica- 
tions, with superconducting Tc above 50 K already been 
achieved [3, 4]. The phase diagram in these materials 
displays both spin density wave (SDW) with the order- 
ing wavevector Q = (tt, 0) and superconducting (SC) or- 
ders [5, 6]. For instance, LaOFeAs is a metal with SDW 
order, which undergoes a first order transition to a su- 
perconductor upon doping [5]. Similarly, CeOi-ajFjjFeAs 
undergoes a doping-tuned continuous transition from the 
SDW to a SC [6]. 

In this Letter, we investigate the interplay between 
these two types of order in terms of symmetry consider- 
ations. It has been suggested that a strong contender for 
the SC order parameter is an "odd-sign" s-wave where 
the SC order parameters on the electron and hole pockets 
acquire opposite signs while there is a full gap at both 
pockets [7-10]. Here, we will show that there is a broad 
class of Hamiltonians with SO (6) symmetry that natu- 
rally give rise to both the "odd-sign" SC order state and 
to an SDW with ordering wave vector Q = (tt, 0). More- 
over, the symmetry analysis reveals a hidden order pa- 
rameter, dubbed "Orbital Density Wave" (ODW). This 
ODW is intimately connected to the SDW and SC orders 
via the SO (6) group; namely they form a six dimensional 
vector under the action of SO (6). The SO (6) symmetry 
is emergent in the sense that it becomes an increasingly 
better symmetry in the limit of longer distances and lower 
energies. We explore experimental consequences of such 
an emergent symmetry in iron pnictides. 

The idea of enhanced symmetry in strongly correlated 
electron systems was first demonstrated in the Hubbard 
model at half-filling [11], and it has been used to describe 
various systems such as the cuprates [12-14] and the or- 
ganic superconductors [15, 16]. The iron pnictides may 
be particularly well-suited for such an approach. The 
Coulomb interaction in these materials may be relatively 
small, as seen in the metallic nature of the magnetic state. 
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FIG. 1: Top panel: Collective excitations in the SC and SDW 
phases. Inside the SC, there are degenerate gapped modes 
Ci.H.z, and a higher energy r] mode. The Q mode, shown 
here as a linearly dispersing Goldstone mode, acquires a gap 
when long-range Coulomb repulsion is taken into account. In 
the SDW state (say, ^i), there are two linearly dispersing 
magnon modes Sy and Sz, and gapped Rx and CX modes. 
The relative size of the and C,x gaps may depend on doping. 
Lower panel: Gap under an applied Zeeman field for different 
collective modes. In the SC, the C,x,y,z modes are split, as 
shown. Inside the SDW phase, the magnon mode Sy obtains 
a gap. All other modes are unaffected to linear order in field. 

This is in contrast to the cuprates, where the interplay 
of SC and magnetism is complicated by Mott physics. 

The emergent S0(6) symmetry leads to several pre- 
dictions for the phase diagram and excitation spectrum. 
The resulting excitation spectra are summarized in Fig.l. 
For instance, inside the SC state, we predict the existence 
of three low energy S = 1 excitations (denoted Q with 
orbital angular momentum L = 0, and charge q = 2 - 
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analogous to the 7r-excitation proposed in the context of 
the cuprates [13] - and a low energy S = excitation 
(denoted 77), which carries L = 1, and q = 2. These 
modes can be probed by inelastic neutron scattering ex- 
periments via 5* = 1 or L = 1 quantum numbers. On 
the other hand, in the SDW state, we predict a low en- 
ergy S = 1 (denoted R) excitation with L = 1, and q = 0, 
which couples to neutrons, and two other excitations (de- 
noted C), with 5 = 1, i = 0, and q = 2, which do not 
couple to neutrons. 

Effective Hamiltonian. — The band structure of the 
iron pnictides involves all five d-orbitals of the Fe atoms. 
Staggering of the out-of plane position of the As atoms 
leads to doubling of the unit cell to include two Fe atoms. 
The band structure has been computed within the local 
density approximation [17] and confirmed experimentally 
in ARPES [18] and quantum oscillation [19] experiments. 
The Fermi surface is composed of two electron pockets 
surrounding the point (0, 0) of the Brillouin zone, and 
two hole pockets around Q = (tt, 0) (or, equivalently, 
(0, tt)). The electron and hole pockets are similar in size, 
leading to an approximate particle-hole symmetry. In 
this analysis we concentrate on one electron pocket and 
one hole pocket, since the physics governing ordering in- 
stabilities is already contained at this level (although the 
other bands are necessary to understand the precise or- 
dering pattern of the SDW and ODW orders within one 
unit cell). 

In a two-pocket model, the most general on-site inter- 
action Hamiltonian is. 
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where, ipjaa creates an electron with spin a at site i. 
a - ■ ■ 6 are pocket indexes: ipj^^ — and = dj^ cre- 
ate electrons near the electron and hole pockets, respec- 
tively. There are five independent coupling constants, 

U4 = U^^, and M5 = U^^; all other interactions are zero. 
For systems with particle- hole symmetry, U4 = U5. The 
Hamiltonian (1) possesses S0(6) symmetry provided that 
the following conditions are satisfied: 



U2 = and U4 



-Ml, 



(2) 



and U3 can take any value. This defines a 2 dimensional 
subspace of S0(6) symmetric Hamiltonians. Remark- 
ably, a one-loop RG analysis of the two band Hubbard 
model with particle-hole symmetry yields an RG that 
flows towards a fixed point satisfying the SO (6) condi- 
tions (2) [8]. The RG is controlled by the size of the 
fermi pockets, and in the limit of infinitesimally small 
pockets the RG flows asymptotically to an SO (6) sym- 
metric point. 

The group S0( 6). — One can show that the interaction 
Hamiltonian at with the constraints (2) commutes with 



the following 15 generators of S0(6), 

Q = -JH(^k,^^k,^ + 4+Q,<,rfk+Q,^-l 



k+Q,(T' 



f}^ = —i 



k+Q,(T'' 



r)' = (f7+)t (3) 



where a takes the values x, y, z, and cr^ are Pauli ma- 
trices. The 15 generators are split into 5 groups: (i) 
the charge operator Q; [ii) three spin operators Sa] {Hi) 
three operators Ra , which measure the difference in spin 
between pockets c and d; [iv) the six operators Cq ; and 
(w) the two operators t)^ . Both Q'^ and 77+ create pairs of 
electrons, [Q,Cq] = ~C^i and [Q,f]+] = -77+, but they 
differ in that t)"*" creates a spin-singlet pair, whereas C,'^ 
creates a spin-triplet pair with zero spin projection along 
the a direction. This difference is reflected in their com- 
mutator with Sa'- [SaTCp] = i^ap-yC,^ and [Sa^f)'^] = 0. 

We can organize the 15 generators (3) by collecting 
them into an antisymmetric 6x6 matrix of operators. 
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whereto = ^{O- +0+) and 3 = ^(O" - 0+), and 
the missing entries are obtained via Lab = —Lha- In the 
matrix form (4) , we can state the commutation relations 
between all 15 operators succinctly: 



cd 



bd 



This is the Lie algebra of the group SO (6). 

S0(6) as a symmetry of the effective Hamiltonian. — 
In order to show that SO (6) is an emergent symmetry of 
the iron pnictides, wc now consider the commutators of 
Lab with the full Hamiltonian, H = Hq + H' , where 



Ho 



E (^^kC^Ck,^ + 4dl+Q.^dk+Q,a) (5) 
kcr 



is the free Hamiltonian. Note that Hq conserves total 
charge Q, total spin Sa, and the relative spin Ra- Hence, 
the only non-trivial commutators of Hq are. 
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In the iron pnictides, the electron and hole pockets have 
similar shapes and sizes [20] . We can then tune the chem- 
ical potential to a point where the system has an approxi- 
mate particlc-holc symmetry, with the low energy disper- 
sion + elj. ~ (nesting condition). At this point the 
free Hamiltonian Hq has an approximate SO (6) symme- 
try. Since H' commutes with the full S0(6) generators. 
H ~ Hq + H' has an approximate SO (6) symmetry. 

In practice, there are three sources of SO (6) symmetry 
breaking. First, the particle and hole pockets do not have 
identical shapes, so that the particle-hole symmetry is 
only approximate. Second, doping the system away from 
the (approximate) particle-hole symmetric point gives. 



(7) 



where jl is the shift of the chemical potential relative to 
the particle-hole symmetric point. In this case, the sym- 
metry of Hq is reduced to S0(4)5flXU(l)Q, the group 
generated by Sa-, Ra, and Q. Increased doping reduces 
the nesting between c and d pockets, and therefore sup- 
presses SDW order in favor of SC. By changing doping, 
we can tune to the S0(6) symmetric point. Finally, there 
are likely to be residual interactions that do not satisfy 
Eq. (2). However, provided that these are not strong, 
the Ci and R modes will remain sharp excitations, and 
the topology of the phase diagram will be unaffected. 

Ground state manifold. — Six different order parame- 
ters are connected by SO (6): superconductivity (A*), 
SDW order (d^,, ^y, and $2), and "Orbital Density 
Wave" (ODW) order (S). For U3 > interactions fa- 
vor the formation of "odd-sign" s-wave superconductiv- 
ity, with a relative minus sign between the electron and 
hole pockets, 

k 

M3 > also favors a SDW instability with momentum Q, 

= J E U.Ad^+Q^^, + 4+Q,.<.'Ck..) . (9) 



Finally, as we will show, it also favors the ODW, 
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k,CT"k+Q,CT - ak+Q,fTCk,a 



(10) 



In real space, S = -^J2ri~'^Y"('^r,c+id-nr^c-id), where 

nr,c±id = 5 J2ai^r-cr ±idra)Hcrcr ±idrcr) ■ ThuS, the ODW 

is a state with alternating c + id and c — id orbitals along 
the X direction - it breaks time reversal symmetry and 
translational symmetry by one lattice spacing along x, 
but preserves the combination of the two. It also pre- 
serves spin rotational symmetry. 



These six order parameters can be organized into a 
vector ha, a — 1 . . .6, where ni = KA, n2 = SA, 713 = 
^x, ^4 = 'J't/i "-5 = ^z, and rig = S. Then, we find 



Lab, fT-c 



(11) 



Hence, the ground state order parameters transform as a 
six dimensional vector under the action of SO (6). 

Phase diagram. — As a first application of S0(6) 
symmetry, we consider the phase diagram of the iron 
pnictides. Symmetry places strong constraints on the 
Ginsburg-Landau (GL) free energy, which in turn con- 
straints the possible topologies of the phase diagram. To 
quartic order in the order parameters, the most general 
GL free energy is. 



Fgi 



d^r 



:(Vn) 



(12) 



Here, the first line is S0(6)-symmetric, while all other 
terms break the symmetry explicitly. The sign and mag- 
nitude of Td can be tuned experimentally through doping. 
The term Fq - not written explicitly - contains all pos- 
sible quartic terms that break SO (6) symmetry, and is 
expected to be small. 

We assume that is small in magnitude, and first 
consider the case rjj = 0. For simplicity, we will also 
assume that F'q is SO(4)fl5 xU(1)q symmetric. A renor- 
malization group analysis of such a free energy was car- 
ried out in Ref. 15. There, the phase diagram was found 
to have three possible topologies, depending on the co- 
efficient g appearing in Fq = 5(|Ap — $2 _ |=;2^2 _|_ 
For 5 > 0, there is a tetracritical point where four phases 
meet: the normal state, SC, magnetism, and a mixed 
state of coexisting magnetism and SC. For <? < 0, there 
is a direct first order transition between SC and mag- 
netism, and first order lines ending at tricritical points 
separating the ordered and normal states. For the spe- 
cial case of full SO (6) symmetry, g = 0, there is a first 
order transition between SC and magnetism ending at 
a bicritical point, with second order lines separating the 
ordered and normal states. 

The precise nature of the magnetic state appearing in 
the phase diagram depends on the sign of r^j. For > Q, 
the magnetic state will be the SDW, whereas for r^j <Q, 
it will be the ODW. Since a ferromagnetic Hund's rule 
coupling gives a positive contribution to r^, we concen- 
trate on the SDW in what follows. 

Low energy excitations. — In order to study the low en- 
ergy excitations that arise due to the approximate SO (6) 
symmetry, we introduce a quantum rotor model [13], 

Hqk = ^ E ^lab - E '^-^^^^ -hY^UM (13) 

i.a<b (ij),a i 
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The model (13) is defined on a coarse-grained lattice, 
such that each site i contains an even number of Fe atoms. 
This is necessary in order to define local operators for 
the SDW and ODW orders. The operators L^ ab and 
n'j satisfy the commutation relations (5) and (11) when- 
ever i = j, and they commute for i =/= j. In addition, 
= for all i and j. In Eq. (13), we allow the 
couplings Ta to be anisotropic, and also include a Zee- 
man field that couples to the spin Sz = Note that 
higher order terms, which are necessary to bound the 
energy from below, are not included explicitly in -ffqR. 
However, in a linear spin- wave analysis such terms only 
lead to weak rcnormalization of the spectra. The collec- 
tive mode spectrum is obtained by solving the Heisenberg 
equations of motion for Li,ab and hf. Figure 1 shows the 
excitation spectrum in the SC and SDW phases. 

Inside the SC, A ^ 0, there are three gapped ( modes 
and a gapped rj mode [21]. In addition, the Q mode, 
shown as a true Goldstone mode of the SC in Fig. 1, is 
charged and acquires a gap of order the plasma frequency 
once Coulomb interactions are taken into account. The 
C modes can be detected in inelastic neutron scattering 
experiments, which probe the dynamical spin structure 



factor xL(k,w) = J^r. 



(5(tJ— i?mo)i where 



(m|5„(k)|0) 

the sum is over all excited states |to) and Emo is the ex- 
cited state energy measured from the ground state. For 
the Cz mode, neutron scattering near the wave vector 
Q yields x"^(k-HQ,w) oc |Ap5(w - (k Q)), where 
u;l{k+q) = |A|2(ri-r5(l-|k|2))/x. The factor of | A 1 2 
reflects the fact that, although is a particle-particle op- 
erator, it can couple to a particle-hole probe such as neu- 
trons inside the SC, where a reservoir of Cooper pairs 
means that particle number is only conserved modulo 
2. Similarly, the rj mode - which connects the SC and 
ODW - can be seen by neutrons with wave vector near 
Q, and the amplitude of the signal is also proportional 
to |Ap. Some of these modes may have been observed 
in a recent inelastic neutron scattering experiment on 
Bao.6Ko.4Fe2As2 [22]. These modes have also been in- 
terpreted as particle-hole bound states in the SC [22]. 
However, unlike these particle-hole bound states, the en- 
ergy of the r] and C modes is independent of temperature, 
a feature that can distinguish between the two proposals. 

Inside the SDW, 7^ 0, there are two gapless magnon 
states, and also (z, Cl, and Rz modes [21]. Unhke the 
superconductor, C| does not couple to neutrons in the 
SDW. The Rz mode, on the other hand, is a particle-hole 
spin-one operator. Therefore, it can naturally be seen in 
spin-polarized neutron scattering. For wave vectors near 
reciprocal lattice vectors G, x'zzO^~^ ^) °^ f{^)^1^{^ ^ 
Lunjk)), where Luljk + Q) = \<^>z\Hr5 ~ re{l ~ \h\^))/x, 
and /(G) is a form factor that vanishes for G = but 
is finite for G 0. Neutrons near Q can also see the Rz 
mode, with scattering amplitude Xzz oc <I'^. 

As shown in the lower panel of Fig. 1, the C modes 



inside the SC are split by an applied Zeeman field, and 
one of the magnon modes in the SDW becomes gapped 
under the field. All other modes arc unaffected to linear 
order in the field. 

Discussion. — In the weak coupling RG analysis of Ref. 
8, the SO (6) fixed point (2) is only reached in the limit of 
infinitesimal pockets. In the iron pnictides, the pockets 
are small but finite. In this case, a functional RG analysis 
finds that the interactions eventually fiow away from the 
SO (6) symmetric fixed point [10]. However, the RG flows 
very close to the S0(6) point, which may dominate the 
finite temperature and energy properties of the system - 
the focus of this work. In some compounds, the size of the 
pockets may be particularly small, and we predict these 
to be the best candidates to observe the SO (6) symmetry. 

In some iron pnictide materials, such as 
LaOi-ajFjjFeAs, the transition between the SDW 
and SC phases is accompanied by a structural transition, 
and is strongly first order [5]. In these materials, 
the SO (6) symmetry may be obscured by the struc- 
tural transition. Fortunately, other materials such as 
Ce0i_2;Fa;FeAs seem to have continuous transitions be- 
tween the SC and SDW phases, without a simultaneous 
structural transition [6]. These materials may provide 
the most promising candidates for the observation of 
S0(6) symmetry. This, in turn, would prove that SC 
and magnetism have a common origin in the electronic 
interactions of the iron pnictides. 

Finally, we discuss the possibility of observing the 
ODW state. Although the ferromagnetic Hund's rule 
coupling favors SDW order over ODW, other residual 
interactions could stabilize the ODW in some of the iron 
pnictides. Elastic neutron scattering experiments can- 
not distinguish between SDW and ODW states, both of 
which give elastic peaks at wave vector Q. On the other 
hand, the ODW does not have low energy magnon states, 
unlike the SDW. Thus, inelastic neutron scattering can 
in principle distinguish between the two states. 
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